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SUMMARY 

A unique method is presented for collecting, studying and interpreting thought-related 
electroencephalogram (EEG) data. The use of a chaos-based nonlinear analysis technique is 
shown to be promising in providing insight into relating conscious thought to specific EEG 
data. A discussion of the practical limitations of this technique is also included. 


INTRODUCTION 

The current primary use of the EEG is in clinical applications. However, because the EEG 
involves conceptually simple instrumentation, is relatively safe, and inherently noninvasive, it is 
well worth investigating for possible use in other applications. Specifically, if it were possible 
for a human operator to reliably control some aspect of his/her EEG, a direct human- 
computer interface could be developed. Thus, the goal of this research is to determine if an 
operator can, with their eyes open, consciously control some EEG feature so as to allow a 
computer to detect one or more bits of information. 

In order to achieve this goal, the laboratory set-up shown in figure 1 was created. The host 
computer is a Macintosh II utilizing National Instruments’ LabVIEW software. The interface 
board is a National Instruments NB-MIO multifunction A/D converter. The EEG amplifier is 
a Grass Model 8-16C EEG Instrument. The electrodes are 9 mm silver-silver chloride contacts 
with collodion adhesive and EEG gel. 

The EEG data to be analyzed were obtained during what is called a "Think Math" test. 
During this test, the subject was instructed to think of a simple math problem (such as 9 x 8 
= 72) whenever the light-emitting diode (LED) was turned on. The periods when the light is 
flashing are called Thought-times. During the off-times of the LED (Gap-times), the subject 
was instructed not to think about mathematics or multiplication. This on/off "Think Math" 
process was repeated for several minutes with the LED flashing periodically two out of every 
eight seconds. The data analyzed in this paper were recorded over what is referred to as the 
"F3 region with respect to Cz" in the 10-20 International System of Electrode Placement 
(figure 2). Data were sampled at 128 samples per second (T = 7.8125 ms), well above the 
Nyquist minimum for data whose highest significant frequency component appears around 15 
Hz. 
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CHAOS ANALYSIS OF EEG DATA 

The application of chaos theoiy to biological systems is not a new concept, nor is the study 
of the chaotic behavior of the EEG. However, the idea of applying chaos theory to determine 
if the EEG can be used as a thought/non-thought discriminator appears to be new. A 
fundamental technique in chaos theory is the reinterpretation of nonlinear time signals as 
multidimensional phase-space geometric objects. The goal is to determine if there exists a 
"strange attractor" in the phase-space diagram of the multidimensional data. A two- 
dimensional phase-space diagram is generated by plotting the original data sequence, x(nT), 
against a time-shifted version of the same data, x(nT + t). Reference 1 suggests that a time 
shift (t) be used which is equal to one-quarter of the period of the most significant frequency 
component in the EEG signal. The reasoning for this value is based on the "mutual 
information content" of a signal. For the analysis in this paper, a time shift of 23.475 
milliseconds was chosen. This corresponds to a dominant frequency of 10.6667 Hz. This is 
intuitively pleasing since the 10.6667 Hz component is in the middle of the alpha region, 
considered by some EEG researchers as the most important EEG frequency region. 


QUANTIFYING THE STRANGE ATTRACTOR 

Figure 3 shows a phase-space diagram for x(nT + r) versus x(nT) for 4 seconds of a 32 
second "Think Math" test. In order to quantify the chaotic strange attractor generated during 
this test, the Standard Deviation/Moment of Inertia (SDMI) transformation was developed. 
This transformation consists of first finding the point-by-point moment of inertia about the 
origin of the phase-space "cloud." This is denoted by 

I x (nT) = x 2 (nT) + x 2 (nT + r). 

The purpose of this transformation is to convert the two-dimensional phase-space data into 
a positive scalar whose weight is a function of the original signal and the time-shifted version. 
After this transformation, a moving-average, standard deviation calculation is performed. This 
is computed as follows: 


(k+1)m 

cr(k) = (1/m) X (I x (nT) - /u. k ) 2 where: 

n=1+mk 


(k+1)m 

li t = x (l/m)I„(nT); 


n=1+mk 


k=0,l,2,...,(N/m)-l 


where N is the total length of the data sequence and m is the standard deviation window length 
(64 points). 

This calculation computes a windowed standard deviation of the moment of inertia of the 
phase-space "cloud" twice per second (every 64 points). Figure 4 shows the results of this 
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analysis by plotting <r(k) and the thought-times versus time. Note that this entire run contains 
four two-second thought times and four six-second gap (non-thought) times. 


CONCLUSIONS 

For this subject, it is clear that the SDMI analysis produces a significant peak in the 
standard deviation plot during the times of transition between thinking and not thinking. In 
attempting to generalize the results shown here, the entire test was repeated on several other 
subjects. These subsequent tests, which utilized either visual thought cueing or audio cueing, 
yielded results very similar to that presented here. Further investigation of the raw EEG did 
indicate the possibility that unwanted artifact (either eye or muscle movement) contributed to 
the peaks generated by the SDMI analysis. Thus, some caution is given in interpreting the 
SDMI peaks as purely thought-related activity. While this research is still in its infancy, it is 
apparent from these results that this technique may be useful for understanding thought- 
related EEG data. 
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Figure 3. Chaotic Phase-Space "Cloud" showing x(nT + r) versus x(nT) for four seconds of 
data. 



Figure 4. Thought-Times and the SDMI EEG data for a 32 second "Think Math" run. 
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